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Do mean-field dynamos in nonrotating turbulent shear-flows exist? 
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Abstract. A plane-shear flow in a fluid with forced turbulence is considered. If the fluid is electrically-conducting then a 
mean electromotive force (EMF) results even without basic rotation and the magnetic diffusivity becomes a highly anisotropic 
tensor. It is checked whether in this case self-excitation of a large-scale magnetic field is possible (so-called W x J-dynamo) 
and the answer is NO. The calculations reveal the cross-stream components of the EMF perpendicular to the mean current 
having the wrong signs, at least for small magnetic Prandtl numbers. After our results numerical simulations with magnetic 
Prandtl number ~ 1 have only a restricted meaning as the Prandtl number dependence of the diffusivity tensor is rather 
strong. 

If, on the other hand, the turbulence field is stratified in the vertical direction then a dynamo-active a-effect is produced. The 
critical magnetic Reynolds number for such a self-excitation in a simple shear flow is slightly above 10 like for the other - 
but much more complicated - flow patterns used in existing dynamo experiments with liquid sodium or gallium. 
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1. Introduction 

One of the most challenging problems in astrophysical fluid 
dynamics is the investigation of the interaction of rotation 
and/or magnetic fields with turbulence. The maintenance of 
differential rotation and the induction of global magnetic 
fields can be the consequences of such interactions. In most of 
the cases, however, the turbulence must be anisotropic and/or 
inhomogeneous to generate interesting global phenomena. 
There is one exception, however, where homogeneous and 
isotropic turbulence interacting with an inhomogeneous mag- 
netic field leads to the well-known fl x J-term in the turbu- 
lent electromotive force (EMF) which together with differen- 
tial rotation can be dynamo-active (Krause & Radler 1980). 
The term, however, vanishes for spectra of the mixing-length 
type (Kitchatinov, Pipin & Rudiger 1994). 

Rotation may not be the only flow whose influence en- 
ables the turbulence to generate global fields. Any vortical 
large-scale motion can be suspected to do the same. The sim- 
plest case to probe the expectation theoretically or in labora- 
tories is, probably, the plane shear flow. 

In the present paper a plane shear flow is considered to 
analyze the main phenomena of the mean-field magnetohy- 



drodynamics and to formulate predictions for its experimen- 
tal realization. We shall show that the interaction of free ho- 
mogeneous and isotropic turbulence with a plane shear flow 
does not lead to the so-called W x J-dynamos (cf. Roga- 
shevskii & Kleeorin 2003). Large-scale dynamos are possible 
only if the turbulence is not uniform along a direction normal 
to the plane of the shear. In this case an a-effect is produced 
by the sheared turbulence which together with the shear itself 
generates global magnetic fields. 

2. EMF of sheared turbulence 

The magnetic-diffusivity tensor rjijk relates the mean electro- 
motive force (EMF) 
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£ = u' x B' (1) 

to the gradients of the mean magnetic field via the relation 
£j = VijkBj.k- (2) 

If the influence of the shear flow is included to first order the 
general structure of the diffusivity tensor is 

Vijk = V0 e ijk + Vl e ijpV'k,p + V2 e i3pUp,k + 

+V3tikpUj, p + riitihpUpj . (3) 

Quasilinear derivations of the coefficients in Q can now be 
performed. The electromotive force Q may be constructed 
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by a perturbation method. The fluctuating fields are repre- 
sented by a series expansion, 

F' = F (0) +F (1) +-F (2) + ... , (4) 

where the upper index shows the order of the contributions in 
terms of the mean shear flow and magnetic field. 

The zero-order terms represent the 'original' isotropic 
turbulence not yet influenced by the shear. The spectral tensor 
for the original turbulence is 



-(o) _ E(k,uj) f kikj 
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(5) 



where the positive-definite spectrum E gives the intensity of 
isotropic fluctuations, i.e. 
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f f E(k,uj)dkduj. 
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We apply the quasilinear approximation (SOCA) to derive 
the higher-order terms in Q- They are generally found by a 
perturbation method from the linearized equations. E.g. the 
linearized momentum equation reads 

» 
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where the upper index shows the order in the mean shear. 
With this equation taken for n — 1, the first-order correction, 
u^\ can be expressed in terms of the given original turbu- 
lence. 

As known, for the isotropic eddy diffusivity 770 one finds 

OOOO 

1 f f r)k 2 E(k,uj) „ , 



(see Riidiger & Hollerbach 2004). The shear-related coeffi- 
cients of (|3j read 
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for i = 1 . . . 4 with the nontrivial kernels 

bv^rfk 12 + uj 2 vr\ (r/ 2 + iv 2 + 61/77) fc 8 
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C3 is negative-definite. Another important property here is 
that the kernels C3 and C4 also exist in the limit v — > 0, i.e. 

8 77 2 fc 4 



C 3 = - 
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For sufficiently small magnetic Prandtl number Pm = v/r) 
the kernel C4 is thus positive-definite. This is not true, how- 
ever, for Pm of the order unity. For Pm = 1 one obtains from 

E3 4 

1 (w 2 -77 2 /c 4 )(5w 2 + 3?7 2 fc 4 ) 
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which has no definite sign. The high-frequency parts of the 
spectrum provide positive contributions to the integral and 
the low-frequency parts provide negative contributions to the 
integrals. White noise (with dE/duj = 0) leads to positive- 
definite values of 774 while the r-approximation (E oc 6(uj)) 
leads to 774 < . Numerical simulations with Pm ~ 1 do 
have thus only a restricted meaning as the Prandtl number 
dependence of 774 is rather strong. 

The so-called r-approximation can be applied as a crude 
representation of nonlinear turbulent effects. The approxima- 
tion replaces the left side of Eq. Q by u\ n /t CO tt, where 
Tcorr = ^coit/u' is the eddy turnover time (£ corl is the cor- 
relation length). The r-approximation can simply be recov- 
ered by substituting uj = and vk 2 = r - ^. In particular, the 
spectrum 

_ (13) 

with l 2 mr /v = r corr can be used to derive well-established 
estimates of the turbulence-induced coefficients l|9} (Vain- 
shtein 1980; Vainshtein & Kitchatinov 1983). 

Only 773 and 774 are important for the simple slab dynamo 
model discussed below. Consider a shear flow with uniform 
vorticity in the vertical z-direction, i.e. 

u y = Sx. (14) 

The shear flow may exist in a turbulence field which does not 
possess any other anisotropy apart from that induced by the 
shear dl4> itself. For experiments in the laboratory it might 
be relevant that the relations 



E(k,uj) = 2u' 2 S(uj)S(k - £-*„), 



£ x = rioDBy 
and 



mSDB x 



(15) 



(16) 



£ y = -r] DB x + T] 3 SDB y 

(D = d/dz) follow from (|2j and Q for a z-dependent field 
imposed in horizontal direction. If the field is imposed in re- 
direction, we have 

e x = -mSDB x , (17) 

and if the field is imposed in the y-direction then 

£ y = mSDB v . (18) 

The sign of 774 is thus opposite to the sign of the expres- 
sion £ X SB XZ and the sign of 773 is the same as the sign of 
£ y SBy :Z . Note that the EMF components are perpendicular 

1 in this case 773 = and no slab dynamo is possible, see Eq. 1301 
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to the mean current J. The standard diffusion-induced EMF 
(without shear) is always parallel or antiparallel to the mean 
current. The EMF due to the standard a-effect is also parallel 
to the mean magnetic field. 

The quasilinear theory provides a positive 774 for small 
Pm and a negative 773 in all cases. We shall see below that 
this constellation does not allow a simple 'W x J' slab dy- 
namo. Our negative result mainly bases on the positive sign 
of 774. Whether 774 is indeed positive must be checked with 
laboratory experiments and/or with numerical simulations for 
various magnetic Prandtl numbers (Brandenburg 2005). 

2.1. Alpha effect 

Now the nondiffusive part of the mean EMF Q is considered. 
If an a-effect exists in the shear flow we have 
S i = a ij B j + .... (19) 

(for more details see Riidiger & Hollerbach 2004). The ten- 
sor a must be a pseudotensor so that an e-tensor has to 
appear in the a-coefficients. The construction of the EMF 
£-i = e ijku'»B' k is the only possibility for the e-tensor to ap- 
pear. Therefore, the subscript of Ei is always a subscript of the 
e-tensor. As the e-tensor is of 3 rd rank an inhomogeneity of 
turbulence with the stratification vector, g = V log u' 2 , must 
also be present for the a-effect to exist. If the shear flow is in- 
cluded to its first order, the general structure of the a-tensor 
is 
ctij = jeijkgk + (otitmu^k + a 2 £ikiuk,j)gi + 

+a^ikig]Ui ik + a^e ik jUi,kgi + a b € i]k Uk,igi- (20) 
If the stratification is along the vertical z-axis it follows from 
<|20} that 

a xx = a 2 g z S = a x S, a yy = -aig z S = a y S, (21) 

and 

19 Z = T. (22) 
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The anisotropy of the a-tensor is described by the difference 
between a x and a y . The so-called turbulence-induced dia- 
magnetism is described by a xy (see Krause & Radler 1980). 
The coefficients of ( I20i read 
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for the pumping term and 
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for two of the kernels. Only the terms occurring in ( 1221 have 
been given. For small Pm, the a x and a y are both of the same 
sign which is opposite to the sign of g z . The a x appears to be 



smaller than the a y . 



3. Slab dynamos 

The shear-flow dynamos without and with a-effect are stud- 
ied with a ID slab model. The dynamo region is finite in the 
z-direction but unlimited in x and 77. The equations for the 
mean magnetic field with zero B z -component read 

8B X 



dt 

dB 



T) D'B X = YDB X - aySDBy - n y SD z B y , 



-^-VoD*B y = 

TDB y + a x SDB x + SB X - r lx SD 2 B x (26) 
with D = d/dz and 

Vx=V4, i h = m- (27) 

The equations can be normalized with 
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dz ' V dz 
d 2 5. 



B X )S 



dz 2 



(29) 



results. The vacuum boundary conditions B X (Q) = B y (Q) 
B x (l) = B y (l) — are applied. 
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Fig. 1. Map for the dynamo instability produced by the W x 
J-effect. Field excitation can be found for sufficiently large 
positive values of the VF-parameter (I30i . 
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3.1. Dynamos without stratification 

Let the preferred direction g be absent in a uniform original 
turbulence. Alpha-effect and pumping of magnetic field both 
vanish in this case. Nevertheless, a hydromagnetic dynamo 
by the W x J-effect might be expected. Figure^shows the 
corresponding stability map. Dynamo excitation indeed ex- 
ists for 

W = fj v (1 + n 2 f, x ) > 0. (30) 

Note that after dlOh r) y = 7/3 vanishes in the r-approximation 
or is otherwise negative-definite. Plane-shear flow dynamos 
are thus not possible in r-approximation - and they are also 
not possible without inclusion of r] x . Only the product r\ x r\ v 
together with the overall shear allows the existence of plane- 
shear flow dynamos. No dynamo linear in the W x J term 
can exist in plane geometry. 

For small Pm, however, r\ x is positive after Eq. Jllh. 
Therefore also for the plane-shear flow (114-i of liquid sodium 
or gallium (Pm < 10~ 5 ) a W x J-slab-dynamo mechanism 
proves to be impossible. It would thus be challenging to probe 
the sign of the r\ x in the laboratory in accordance with the 
Eq. {T7J, 

Another situation holds for the case of magnetic Prandtl 
number of order unity which is used in most of the numeri- 
cal MHD simulations. As outlined above, the 7/4 can then be 
negative if the frequency spectrum of the turbulence is steep 
enough (if it is too steep then 773 vanishes!). A limiting mag- 
netic Prandtl number between the two cases must exist and 
can be fixed by numerical integrations. 



time for the case of isotropic a-effect. In both considered 
cases a typical condition for dynamo excitation is 



a y S> 10. 



(31) 



(Fig. 0. A reformulation of this relation with S ~ u/L yields 



Rm > 10- 



L 1 



a y H 



(32) 



where Rm=L7j/?7o is the magnetic Reynolds number of the 
mean flow defined with the eddy diffusivity (jSJl. The estimates 
770 ~ u'^corr/3 and ay ~ £^ OTI /H then reduce the inequality 
d32t to the condition 



U < 



31 



(33) 



The fluctuating velocity should be sufficiently small to pro- 
duce a dynamo. It cannot be too small, however, as the 
eddy diffusivity must remain larger than the microscopic one, 
770 > 77. This finally yields the condition 



V 



< u < 



U£ r 



31 



(34) 



which must be fulfilled for dynamo self-excitation. The mi- 
croscopic magnetic diffusivity 77 of liquid sodium or gallium 
is about 0.1 m 2 /s. If the channel width, L, is as large as 1 m, 
then for the most optimistic case (£ corl ~ L) the range (I34> 
can be realized with shear velocities exceeding 1 m/s which 
is of the same order as in other dynamo experiments with a 
much more complicated geometry. 
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Fig. 2. Stability map for a global dynamo solution of 
Eqs. d29t for T — —0.5. Magnetic fields are excited in the 
parameters region above the lines. Solid: a x = a y , dashed: 
a x = 0.1d y . 



3.2. Dynamos with stratification 

The a-effect dynamos may be considered now ignoring the 
anisotropic dissipation in the Eqs. j29> . For given shear, self- 
excitation of magnetic fields can be observed for sufficiently 
large a y and/or large enough shear. All the solutions are os- 
cillatory. The oscillation time is of the order of the diffusion 



4. Summary 

It is suggested that the basic effects of the theory of the turbu- 
lent dynamo which are usually concerned as special proper- 
ties of rotating fluids can also be found for the plane shear 
flow. The same is probably true for any flow with global 
vorticity. The elementary structure of the shear flow largely 
simplifies the consideration. In the Appendix we have shown 
with the same concept as above that the generation of large- 
scale vorticity by a uniformly sheared turbulence is not pos- 
sible and the generation of magnetic fields by such a flow can 
only be hoped for large magnetic Prandtl numbers. 

The dynamo instability can be realized, however, when 
the turbulence is not uniform. Then the stratified turbulence 
produces the a-effect which can excite an oscillatory mean 
magnetic field in the shear flow. This opens a possibility for 
the realization of a turbulent dynamo in the laboratory with a 
quite simple flow geometry - if a nonuniform turbulence can 
be created in the channel flow. Estimates of the excitation 
condition for such a dynamo, however, shows it to be at the 
limit of current experimental possibilities. 
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Appendix A: Hydrodynamic stability of the 
shear flow 

It is important for the above consideration that the linear 
shear flow (1141 is hydrodynamically stable under the pres- 
ence of nonlinear shear terms in the correlation tensor Qij . 
This question has been addressed by Elperin, Kleeorin and 
Rogashevskii (2003). With a dispersion relation formulated 
on the basis of JA4I below) an instability has been constructed 
for plane wave disturbances with spatial inhomogeneity in 
the vertical z-direction. In contrast to this, we shall show that 
within the first-order smoothing approximation the shear flow 
is stable. 

The one-point correlation tensor 

(Al) 



Q tj = u'i(x,t)u'j(x,t) 
in its linear form reads 

Qij = P$ij ~ v (U%,j 



*J,V 



(A2) 



Here vq is the isotropic eddy viscosity; the turbulent pressure, 
P, includes all coefficients of the Kronecker tensor 8ij . 

The experiment by Champagne et al. (1970) with sheared 
turbulence indicates that the linear relation (|A2ji cannot be the 
whole truth. In the experiment the rms downstream velocities 
are systematically larger than the rms velocities in the cross- 
stream direction. The turbulence intensities for these two di- 
rections should, however, be equal after Eq. (|A2ji which can 
be read as 



,/2 _ 



P, u 
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P, u' x u' = —v S. 
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The same remains true if higher-order derivatives such as 
Uijn + Uj.m are included. 

However, if the mean shear is indeed the only reason for 
anisotropy, one has also to involve nonlinear terms, i.e. 

Qij = PSij - i/ (itij + %i) + v\u ltk u^ k + 

+P2Uk,iUk,j + vz (u it kUk,j + %fcUfc,i) ■ (A4) 

By this expression the horizontal intensities can differ if the 
coefficients v\ and vi of the nonlinear terms do not coincide, 



C - K 2 = {yx - "2)s 2 
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It should be v\ > V2 for agreement with the aforementioned 
shear-flow experiment. 



The first-order term of Eq. JA2> is the eddy viscosity 
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(Krause & Rudiger 1974). The second-order correction to the 
correlation tensor reproduce the nonlinear terms of JA4> . The 
coefficients v\, v-i and v% result as 



Ki(k, u))E(k, uj)dk dto, n = 1 . . . 3 , 
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with the kernels 
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105(w 2 + ^ 2 fc 4 f 
All the kernels have negative (third) terms in their numer- 
ators. Nevertheless, all the coefficients JA7I > are 'almost al- 
ways' positive. They are positive definite in the most popular 
simplifying cases, i.e. 

- within the r-approximation: 

5 ^vr 32 



v\ = v-i = 



21 



,(0)2, 



U-2 



are always positive 

- for white-noise spectra: 
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dk 

' = Wv J E ^¥ (A10) 
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are always positive. The gap between \A9\ and ( lAlOt is filled 

by 

2 w 
E(k,u) = j— —E(k). (All) 

The coefficients (IA7> are positive-definite with (IA1H . 

For small disturbances, u, of the mean flow depending on 
z the cross correlations defined by \AA\ read 
Qxz = -v Q Du x + v 2 SDu yi 

Qyz = -v Du y + v 3 SDu x . (A12) 

The linear equation system for the disturbances is 

—^ - v D 2 u x + i/ 2 SD u y = 0, 

di 



— ' Su x — vqD Uy + v^SD u x = 0. 



(A13) 



dt 
It reduces to 

{v2V Z S 2 - vl) D 2 u x + v 2 S 2 u x = (A14) 

in the stationary case. The ID problem should be accom- 
plished by boundary conditions imposed at (say) z = and 
z = H. But Eq. JA14t possesses a solution only if 



H = n\ 



'vl-V2V 3 S 2 



-V2S 2 
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irrespectively of whether no-slip or stress-free boundary con- 
ditions are applied. The second term in the numerator of 
( IA15J must be small compared to the first one. Otherwise 
the nonlinear correlations JA4> which neglect the third and 
higher order terms in the mean shear cannot be applied. 
Therefore, an instability only exists for i>2 < 0. The quantity 
V2 proved to be positive for almost all spectra. The shear flow 
dl4> thus proves to be stable in the hydrodynamic regime. 



